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The vicinity of the critical point of the three-state Potts model on a Kagome lattice is studied 
by mean of Random Matrix Theory. Strong evidence that the critical point is integrable is given. 



I. INTRODUCTION 

The q-state Potts model Q, being one of the simplest generalization of the Ising model, is a challenging model 
of classical lattice statistical mechanics. Despite its apparent simplicity it has not yet been analytically solved for a 
number of states q > 2 even in two dimensions (for a review see Q]). The square, the triangular and the hexagonal 
lattice have been shown to be integrable at the critical point. The free energy the correlation function have been 
calculated. Out of the critical point the free energy has not been calculated, moreover strong indications that the 
model is not integrable have been given On another hand the two dimensional Kagome lattice is of great 
interest both from the statistical point of view and from solid state point of view. The anti ferromagnetic Kagome 
lattice provides an example of geometrically frustrated spin system M. Several experiments on the so-called SCGO 
compounds family suggest the existence of a " spin-liquid" phase transition Q . It has also been seen as topological 
spin glass M. 

The critical point of the Potts model on the Kagome lattice is known in the Ising (q = 2) case only. However, an 
algebraic variety has been conjectured to be critical for the checkerboard square lattice with multi spin interaction in 
JjJ . From this conjecture the critical point of the Potts Kagome lattice has been deduced . This value has recently 
been confirmed with high precision by extensive Monte-Carlo simulation so that the actual critical point is at 
least extremely closed from the conjectured point. The aim of this letter is not give another numerical confirmation 
of this conjecture, but to address the question of a possible integrability of the critical point. Indeed it is possible 
that the critical point of the Kagome q = 3 Potts model is of the same nature as the critical point of the square, 
hexagonal or triangular lattice, ie an integrable point. It is however also possible that this point is not integrable, as 
is, for example, the the critical point of the three dimensional Ising cubic lattice This letter is intended to give 
an answer to this question. It is organized as follows: in a first paragraph we sketch briefly the basics of Random 
Matrix Theory (RMT) applied to classical lattice statistical mechanics, and recall how it acts as an 'integrability 
detector'. In the next paragraph we give a more convenient formulation of the Kagome lattice as a checkerboard 
Interaction-Round-a-Face (IRF) model, and we discuss the symmetries of the corresponding transfer matrix. In the 
final paragraph the results are presented and discussed. 

II. THE RMT METHODOLOGY 

The RMT lattice statistical mechanics has first been introduced for Quantum one dimensional system P|lC|] , it has 
then been extended to classical two dimensional lattices jllj| . Roughly speaking the main idea is to compare the given 
Hamiltonian to an 'average Hamiltonian'. Among all Hamiltonians, integrable Hamiltonians are extremely peculiar 
and show significant difference with 'average Hamiltonian'. To be more specific one is interested in the statistical 
properties of the discrete spectrum viewed as an infinite set of real numbers. One compares the statistical properties of 
the spectrum of the system under consideration with the one of random matrices from a suitable statistical ensemble. 
For an integrable system it is found that the spectrum as many properties of a set of independent numbers, ie the 
suitable matrix statistical ensemble is the set of the diagonal matrices with independent normal centered entries. By 
contrast, for a non- integrable time reversal symmetric system the spectrum is well described by spectrum of matrices 
from the so-called Gaussian Orthogonal Ensemble (GOE) of symmetric matrices with independent normal centered 
entries. The eigenvalue spacing distribution P(s) discriminates well between these two limiting cases: P(s) ds is the 
probability that the difference between two consecutive eigenvalues belongs to the interval [s, s + ds]. It is simple to 
show that for independent eigenvalues one has P(s) = e~ s . For GOE matrices P(s) is extremely close (actually strictly 
equal for 2 x 2 matrices) to the Wigner surmise W(s) = -|s exp(— 7rs 2 /4). Provided that the proper normalization 
explained below has been performed, and if the system size is large enough, then the distribution P(s) of the system 
under consideration is close to one of the exponential or Wigner distribution. To quantify this "proximity" one 
introduces the parametrized distribution 
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which is one of the infinitely many possible interpolations between the exponential and 
Wigncr distribution. For /3 = one recovers the exponential law and for = 1 one recovers the Wigner law. The full 
lines on Fig. |^ represents these two distributions. Parameter is interpreted as the repulsion between eigenvalues. 
P(s) is a quantity involving only pairs of eigenvalues; a quantity involving more than two eigenvalues is the spectral 
rigidity 
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A 3 (A) = ( - min / (N u (e) - ae - bfde ) , (2) 
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where N u (e) = 0(e — is the integrated density of unfolded eigenvalues and (. . .) Q denotes an average over a. It 
provides a finer test of the closeness of a spectrum with a GOE spectrum or with independent numbers. These two 
limiting cases are shown of Fig. || as full lines. For the classical lattice system, like spin models or vertex models, the 
transfer matrix [ [l3[ is considered instead of the Hamiltonian which generally has a trivial spectrum. 

The operational procedure has been explained in details in Jl2| , p~3| | . It consists in two main steps. One first sort 
the states according to their symmetries: states of different symmetries have to be considered separately. Then, since 
universal properties cannot be found on the row spectrum, a procedure of "normalization" , known as the unfolding, is 
performed: it consists in imposing a local density of eigenvalues of the order of unity. Then the statistical properties 
of the sorted and unfolded spectrum are analyzed. Note that this analysis requires real eigenvalues. This is generally 
insured by the of the Hamiltonian. 



III. TRANSFER MATRIX OF THE POTTS MODEL ON THE KAGOME LATTICE 

The Kagome lattice is a two dimensional lattice depicted on Figure [l] The black circles represent the sites, and the 
lines are the coupling bonds. The partition function reads 
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< i,j > denotes the neighboring sites on the Kagome lattice, 8 is the Kronecker symbol, Kij = K\,Ki or K3 for 
each of the three types of bonds, and Oi can take one of the three values 0, 1, 2. The Kagome lattice can be regarded 
as a lattice of basic "cells" (1,2,3,4) as shown on Fig. [I]. Each cell consists in a up pointing triangle sharing a site 
(site labelled on the figure) with a down pointing triangle. These cells are disposed on every other (say black) 
plaqucttcs of checkerboard. To see the correspondence with an IRF model on a checkerboard lattice, one perform a 
partial summation on the central site (site 0), seeking coupling constants such that 

W / IR F (CT1, (72, 03, CT 4 ) = ^ W/ kag(o-l, 0"2, ^3, CT 4 ; CT ) (4) 

This is possible if one introduces multi spin interaction. The general correspondence is given in Appendix |a|. For the 
isotropic K\ — K2 = K3 = K case the Boltzmann weights of the plaquettes are given by 

W ka . s (a 1 ,a 2 ,(T 3 , a 4 ; a ) = K s °^ +s ^ +s ^ +5 °<*+ 5 ^+ s ^ (5) 

and 

WW(cti, (T2, 0-3, 0-4) = L S ^ 2+&3A L 5 v 1 - 3+S2A L S d 1A+52 ' 3 T s ^> 3+S ^>* +S ^>* +S2 ^ F Sl - 2 - 3A (6) 

where 5i t j = 8(jJi,crj) and 8i_j t k = 8(<Ti, <jj , o^) (cj = 1,2,3,4). The parameters of the IRF model are given as a 
function of the parameters of the Kagome by: 

L h = KL V = KL d = 

(K + 2)(K 3 + K+1)K 

(2K 2 + 1) 2 {) 
_ (K 4 + 2){2K 2 + K + l) 3 
F ~ (K 3 + K + l) 4 
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FIG. 1. Correspondence between the Kagome lattice and the checkerboard IRF model. 



From a computational point of view it is more convenient to use the IRF representation, this is the form chosen in 
the calculations presented here. A first difficulty in constructing the transfer matrix comes from the fact the actual 
transfer matrix which brings from row A to row C on Fig. [I] is the product of the two transfer matrices bringing 
respectively from row A to row B and from row B to row C. As a result the calculation of a single entry of the transfer 
matrix involves a combination of sums and products over N terms, N = 3 L being the size of the matrix (L is the 
linear size of the rows, L = 6 on Fig. This greatly slows down the calculations. Another difficulty comes with the 
fact that, unlike simpler lattice, there is no way to make symmetric the transfer matrix T, so that the eigenvalues 
are complex. It is still possible to deal with complex eigenvalues. In this case integrability also leads to independent 
eigenvalues, and non integrability leads to repulsion. However, in practice, it is much more difficult to make the 
distinction, since in two dimensions repulsion is less effective. Therefore we have chosen to analyze the statistical 
properties of the operator T = (T + T)/2. 

In the actual computation, the transfer matrix calculated is the diagonal-to-diagonal transfer matrix. As explained 
in the previous paragraph one needs to find the symmetries of the operator T. Let us define the translation operator 
i, the reflection r, and the color operator c 

t Wqi&Ii • ■ ■ >0X-1 > = |fl, cr 2 , • • • (To > (8) 

r \a ,ai,- ■ ■ ,a L -i > = \a L -i, ctl-2, • • • , Co > (9) 
c \a ,ax, • ■ • , > = \a + 1, cri + 1, • • • , ol-i + f > (fO) 

where the sum are taken modulo 3. It is easy to see the three operators t, r and c commute with T, and that c 
commutes with both t and r. Thus the symmetry group is the direct product T>l ® T>3 where T>l is the dihedral 
group of index L. The irreducible representations of this group are simply the matrix product of elements of T>l and 



2?3 (see details in |fj|). The sizes of the blocks corresponding to the different representations are given in Tab. IV 
In the statistical analysis the blocks of size smaller than 100 have been discarded to minimize small size effects. 



IV. RESULTS AND CONCLUSION 



Using the IRF representation of the Kagome lattice we have calculated for different values of the exponential of the 
coupling constant K the symmetrized diagonal-to-diagonal transfer matrix T{K) for L — 8. The size of the matrix 
T is 3 s = 6561. For sizes smaller than L = 8, the blocks are quite small, giving poor statistics. On another hand the 
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size L = 9 is out of reach of your computing possibility, especially due to the facts that T is a matrix product with all 
entries non zero (in contrast with most quantum Hamiltonian) as explained in the previous paragraph. So we stick to 
the case L = 8 and perform the analysis for different values of K. The matrix is projected in the different invariant 
subspaces, yielding matrix blocks of smaller sizes (see Tab. 0). Within blocks there is no degeneracy left, proving that 
the symmetry group used is the largest symmetry group for this model, as it should. The spectrum of the 13 larger 
blocks is then unfolded, leading to 2187 eigenvalues distributed in 13 sub-sets. The eigenvalues spacing distribution 
P{s) as well as the rigidity A3 are calculated. We present on Fig. || and Fig. || the results for two values of K. Firstly 
for K = K c where K c is the value predicted by the Wu's conjecture, which is known to be at least extremely close 
the the actual critical value. And secondly for K = 2 which is sufficiently far from the critical value, and neither too 
small nor too big: extremal values leads to accuracy problems. The results are unambiguous: close to the transition 
point the distribution of eigenvalues is very close to an exponential, whereas out of this point it is close to a Wigner 
distribution. This is clearly seen on Fig. g where the numerical data are presented for both cases together with the 
exponential law and the Wigner distribution. The best fit values of [3 (see Eq. [l]) are /3(K C ) — 0.01 and /3(2) = 1.02. 
This is a strong indication that the critical point is integrable. To further support this result, the rigidity Eq. || has 
also been calculated for the same value of K. The results are shown on Fig. ||, where the numerical data and the two 
limiting cases of a GOE spectrum and of independent numbers are presented. The agreement between the expected 
behavior and the observed behavior is very good, at least up to an "energy scale" of the order of 15, (ie 15 eigenvalues 

since the mean density of eigenvalues is one). The number variance £ 2 (Z) = ([N u (e + 1/2) — N u (e — 1/2) — Z] 2 ^} 

where the brackets denote an averaging over e has also been calculated, giving results not presented here, in perfect 
agreement with the results for the spectral rigidity A3. 

In conclusion a RMT analysis performed on the symmetrized transfer matrix of the Kagome lattice seen as an IRF 
checkerboard model shows that the critical point, very close or equal to the Wu's conjecture, is an integrable point. 
This is in agreement the conformal theories of the two dimensional lattice statistical mechanics. The only example 
of a critical point where the transfer matrix eigenvalues spacing distribution is of Wigner type is the 3-dimcnsional 
Ising case. 
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TABLE I. Sizes and degeneracies of the different blocks for L = 8. The first row is the label of the representation, the second 
is the dimension of the representation (ie degeneracy of the corresponding eigenvalues) and the third row is the size of the block 
to be diagonalized. 



4 



~£2- 

Q_ 



0.8 



0.6 



0.4 



0.2 




FIG. 2. Distribution of eigenvalue spacings for K ~ K c . The linear size of the lattice is L = 8, leading to Hilbert space 
of size 6561. 1910 eigenvalues are included in the statistic. The diamonds corresponds to the critical point K = K c , and the 
square corresponds to the the value K — 2. 




FIG. 3. Rigidity A3 as a function of A for the same data as in Fig. [| 
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APPENDIX A: KAGOME LATTICE AND CHECKERBOARD IRF MODEL 

For a non-isotropic Kagome lattice the Boltzmann weight of a basic cell (see Fig [l]) is 

H / kag(CTi, <J 2 , cr 3 , 04, cr ) = K 1 K 2 K 3 (Al) 

and for IRF checkerboard it is 

ir/ / \ r 5l. 2 Tfe, 4 r^l, 3 r*2, 4 7-^1,4 7-^2,3 7-^1,2,3 r ^1.2,4 r^l, 3,4 rfe, 3, 4 r^l,2,3,4 / r 34\*12^34 / r 24\^13 , '24 / A r>\ 

H / IRF(CTl,CT2,Cr 3 ,(T4) — L 12 L 34 L l3 L 2 ± h u L 23 L l23 L l2i L 134 L 234 -L 1234 (^i 2 J \ L 13) \ AZ ) 

The notation are the same as in the text Eq. ||. The parameters of the IRF model are given as a function of the 
parameters of the Kagome by: 



L l2 {K x ,K 2 , K 3 ) = L M (K l7 K 2 ,K 3 ) = K x L l3 (K 2 ,K z ) = K X L 2A {K 2 ,K 3 ) 



(K 2 K 3 + K 2 + K 3 ){1 + Kl + Kl) 



2Kl + K^K^ + AK 2 K 3 + 2K% 
L U {K 2 ,K 3 )=L 23 (K 3 ,K 2 ) 



G 



K 2 2 +Kj + 1 
2K 2 + 2Kl 

L 234 (K 2 , K 3 ) = L 134 (K 3 , K 2 ) = L 124 (K 3 , K 2 ) = L 123 (K 2 ,K 3 ) 
= (Kj + 2K 2 )(K 2 Kj +K 2 + \){Kl + 2K 3 f 
{{K 2 K 3 + K 2 + K 3 )(Kl + K\ + 1))2 
L\\{K 2 ,K 3 ) = L\\{K 2 ,K 3 ) 

_ (2K 2 K 3 + l)(2Ki + K\Kl + AK 2 K 3 + 2Jf|) 



L 1234 {K2,K3) = 



{K\ + K\ + 1)(X 2 ^ 3 + K 2 + X 3 ) 2 
(K 2 2 K 3 2 + 2)(X 2 2 + g 3 2 + l) 4 
((if 2 + 2if 3 )(^3 + 2X 2 )) 2 

(if 2 if 3 + jf 3 + if 2 ) 4 

X ((K$K 3 + K 3 + l)(K 2 Kl + K 2 + l)) 2 
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